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The low-temperature generalization of the mode-coupling equations corresponds to the 
dynamics of mean-field disordered models in the glassy phase. The system never achieves 
equilibrium, preserving the memory of the time elapsed after the quench throughout its 
evolution. 

A concept of effective temperature can be made quite rigorous in this context by con- 
sidering readings of thermometers in different time-scales and the thermalization of weakly 
coupled subsystems. 

In the past few years progress has been made in the analytical understanding 
of glassy dynamics. On the one hand, it was realized that the out of equilibrium 
dynamics of mean-field disordered models was solvable for long times, and that the 
hition showed aging phenomena El qualitatively close to the ones of real systems 



>@). 

On the other hand, it became clear that, at least at the mean-field level, the 
important distinction is not between models with and without quenched random- 
ness cf'Ef but rather between different kinds of nonequilibrium dynamics. This be- 
came particularly evident when it was noticed that the mode-coupling equations 
for supercooled liquids EP could be obtained as the dynamical equations for mean- 
field disordered models in the high-temperature .phase 0). By considering the low- 
temperature dynamics of these spin-glass modelsliP, a low-temperature extension of 
the mode-coupling equations was immediately obtained 

Indeed, except for the crucial (for structural glasses) question of the behaviour 
near the glass transition, and the related problem of cooling-rate dependence, the 
models described above yield a rather realistic picture of glassy dynamics well below 

Let us consider n modes 0±, O n , their correlations C a b(t, t w ) = {O a (t)Ob(t w )) 
and their responses to perturbations h a conjugate to O a : 



SO a (t) 



, Xab(t,t w )= / dt' R ab (t,t') . (1) 

h=0 JU> 



5h h {t v 

One can the write, in general, Schwinger-Dyson equations for their correlations 
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C = (C a b) and responses R = (R, 
dC ab (t 

1 L w, 



Ot 



= 2TR ab (t w ,t) + (-Vac{t)C cb {t, t w ) + jf " dt"D ac (t, t")R cb (t w ,t' ) 

+ f dt"U a c(t,t")C cb (t",t w ) , (2) 
- Jo 

0B "'" (L / '" ) - S(t - t w )8 ab + £ (-/U^(M™) + T dfEUt^Rcbit",^ 



(3) 



In mean-field models one can close the Schwinger-Dyson equations into a set of 
dynamical equations involving only the correlation and response functions: 

D ab (t, = F ab (C(t, t')) , E ab (t, t') = Fab,cd{ C ^ t')) R cd , (4) 

c,d 

dF d 2 F 
Fabil) = t— , F ab , cd {q) = , (5) 

(J (lab VqabOqcd 

where F is a given function determined by the model. & 

If the system equilibrates, 1) the two-time functions become time-translational 
invariant (TTI), 2) there is reciprocity C ab (t — t w ) = C ba (t — t w ), and 3) we have the 
fluctuation-dissipation theorem (FDT): 

R ab {t - t w ) = -^(t - t w ) , xu(t ~t w ) = - (Cia(O) - C 12 (t - t w )) . (6) 
T ot w T 

Putting this information in Eqs. (2) and (3) one obtains a single equation for the 
correlations: 



dC ab (t-t w ) ^ 1 



oc 1 



^ /^ac C c fe(t — t w ) + — [-Dac(O) C cb (t w - t w ) - D^.C cb 



OT „ T 



+^f t dfD ac (t-f) dCcb ^' tw) , (7) 



ot" 

where D™,C%£ stand for lim^oo C ac (t) and D^. = lim^oo EL c (t), respectively. 
Equation (7) is the usual equilibrium MCT equation for liquids. B It would be quite 
difficult to guess (2) and (3) from (7)! 

As an example, consider the single-mode case 

F(q) = <f . (8) 

The mass fi(t) is chosen so as to impose normalization at equal times of the auto- 
correlation C(t,t) = 1. With these choices Eq. (7) is the simplest Mode^oupling 
equation proposed by Leutheusser and Bengtzelius, Gotze and Sjolander.liirEP. The 
two-time dynamical equations correspond to the low temperature dynamics of a 
spin- glass model introduced by Crisanti and Sommers.© 
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Fig. 1. Decay of the auto-correlation function in the high temperature phase. The plateau develops 



The correlation decay in the high-temperature phase is obtained by solving (7) 
where one sees the familiar a and (3 decay (see Fig. 1). 

As one approches the transition temperature from above the plateau in qEA 
becomes larger and larger, but at the same time something else happens: going back 
to the original out of equilibrium equations (2) and (3) and solving them starting 
from a random initial condition we discover that the time needed to equilibrate 
diverges at the transition. Hence, for temperatures at or below T g we no longer can 
assume TTI or FDT, and Eq. (7) becomes irrelevant. 

In fact, the two-step process in the correlation decay becomes a waiting-time 
dependent two-step process: the correlations fall below qEA ever more slowly as one 
considers an older system (see Fig. 2). Non-trivial T-dependent exponentsEif a and 
(3 characterize the relaxation around qEA- 

We thus see that one of the equilibrium conditions (TTI) is violated: we might 
expect that also the FDT will be violated. Indeed, this is so: a parametric plot 
of x(t,t w ) vs. C(t,t w ) does not yield a straight line with gradient — 1/T as in 
equilibrium (cfr. Eq. @), but a family of curves as in Fig. 3. As one considers larger 
t w the curves approach two straight lines©, one with gradient —1/T, corresponding 
to large and similar times and another with gradient —X/T (X < 1), corresponding 
to large and widely separated times (i.e. the aging regime). 

The fact that T/X — the fluctuation-dissipation ratio — might be related to an 
effective temperature T q q = T/X was noted several times, in particular by Hohen- 
berg and Shraiman in the context of spatiotemporal chaos and weak turbulence. 

Recently M> , it was argued that indeed T e g deserves the name temperature in 
that 1) it is related to the reading of a thermometer that is brought in contact 



at C — qsA 
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with the glass, 2) it decides the direction of heat-flow and 3) it is a criterion for 
equilibration. 

Let us first consider a very simple thermometer consisting of a harmonic oscil- 
lator that is weakly coupled to an observable O of the glass. For example, in Fig. 4 
we show how this can be done if O is the magnetization. 

In the absence of coupling 0(t) has fluctuations with, we assume, zero mean 
(0(t)) = and correlations (0(t)0(t w )) = C (t,t w ) = O(N). The response of the 
system to a field h conjugate to O is Ro(t,t w ) = 5{0(t)) /5h(t w )\h=o- 

The oscillator takes up energy from the fluctuations of O, and dissipates it 
through the response of the system until a stationary regime is achieved. If the 
system is in equilibrium, equipartition of energy implies that T = (E osc ) (we have 
set Boltzmann constant to one). If we now take To(u!,t w ) = (E osc ) as the natural 
definition of frequency and time dependent temperature of O, a simple calculation 
yields: 

T (u,t w )^(E 0SC ) tw = ^^ , (9) 

Xoy^o, t w ) 

where we have used the waiting-time dependent susceptibility and correlations de- 
fined from: 

\x' (u>,t) + ix" {oo,t)] exp(iu;t) = / dt' R(t, t') exp(zwt') , 

Jo 

C'(uj,t) + iC"(uj,t)\ exp(iut)= [ dt' C(t,t') exp(iwt') . (10) 
J Jo 

If a system is in equilibrium, FDT holds and To(oo, t w ) is independent of the ob- 
servable O, the waiting time t w and the frequency to, and is equal to the temperature 
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Fig. 3. The susceptibility xi^^w) vs - tne auto-correlation function C(t,t w ) at T < T g . The full 
curves correspond to different total times t, from bottom to top, t = 12.5, 25, 37.5, 50, 75. The 
dots represent the analytical solution when t w — > oo. 
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Fig. 4. An effective temperature measurement for a magnetic system. The coil is wound around 
the sample, which is in contact with the bath. Th 'thermometer' is the L — C circuit. 
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of the bath. 

In the system defined above, we obtain for large frequencies To(oJ,t w ) = T 
as uj — > oo. If instead we consider a frequency- waiting time domain such that 
C(t w + ^,t w ) < he A (i-e. we are probing the aging scale), then: lb (a;, t w ) = T/X > 
T. 

Consider now an experiment in which we connect the oscillator to an observable 
0\ and let it equilibrate at the temperature To 1 (uJ,t w ), after that we disconnect it 
and connect it to another observable 2 and let it equilibrate at the temperature 
To 2 (w,t w ). The net result is that an amount of energy To 1 (u>,t w ) — To 2 (^>,t w ) was 
transferred from the degrees of freedom associated with 0\ to those associated with 
2 : the flow goes from high to low temperatures. 

This is somewhat like touching two points of a glass that has been thermalizing 
for a long time with a copper wire. Even if the glass is still out of equilibrium we 
would be surprised if one could obtain an energy flow through the wire this way. 
Proposing that this cannot happen is equivalent to saying that for an 'old' glass, 
different observables should have in the same frequency range the same temperature. 

In order to test this idea, we enlarge the model we have been discussing by 
considering two modes, with F given by: 

F(q) = q p u + KV22 • (11) 

We impose normalization at equal times of the autocorrelation of both modes: 

C u (t,t) = C 22 (t,t) = l. (12) 

through two Lagrange multipliers fJ,n(t), fi 2 2 (i)- We couple both modes through 
M12 = M21 = e (cfr. Eqs. (2) and (3)). The constant K ~ 0.7 was introduced in order 
to break the symmetry between the modes. 

Remarkably, it turns out that one can close the equations with two ansatze for 
the long-time aging behaviour. In terms of the effective temperatures their meaning 
is: 

1. Thermalized aging regime. The effective temperatures associated with the ob- 
servables 0\, 2 are equal to each other for frequencies and waiting times in 
the aging regime but different from the temperature of the bath. At higher 
frequencies, their temperatures coincide with the one of the bath: 

Ti(oj, t w ) = T 2 (lo, t w )=T , t w -> 00 , C a b> Qab A > quasieq. , 

Ti(w, t w ) = T 2 (uj, t w ) / T , u -> , t w -> 00 , C ab < q^ b A , aging . 

Not surprisingly, in this case we find that 0\ and 2 are strongly coupled (also 
in the aging regime) in the sense that the cross responses in the aging regime 

Ri2(t,t w ) = ' Rn{t,t w ) = Y~Qj~ ' ( 13 ) 



are of the the same order of the self response functions (X > 0). 
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2. Unthermalized aging regime. The effective temperatures associated with the 
observables Oi, 2 for combinations of frequencies and waiting times corre- 
sponding to the aging regime are neither equal to each other nor to that of the 
bath, while for higher frequencies they both coincide with the one of the bath. 



Ti(u,t w ) = T 2 (uj,t w ) = T , t w — > 00 , 
Ti(u, t w ) / T 2 (lo, t w ) , u^0,t w 



00 , 



C a b > C, 



Cab < C^ A 



EA 
ab 
■EA 
ab 



quasieq. 
aging . 



In this case, 0\ and 2 are effectively uncoupled (in the aging regime): 



R\2{t, t u 



X12 dC\2 

t dt„ 



R2l(t, t w ) 



X21 dC2i 
T dL„ 



(14) 



with X12 — > and X21 — > 0. 

In Fig. 5 we plot x(C) for the two uncoupled systems (e = 0) evolving from the 
initial condition C u (0, 0) = C 22 (0, 0) = 1 and C 12 (0, 0) = C 2 i(0, 0) = 0. We see that 
Xn ^ X22 while X12 = X21 = (unthermalized case). In Fig. 6 we consider the 
same two systems, this time coupled weakly (e = 0.7) Clearly, after a short transient 
associated to short times, all curves Xij(Cij) become parallel. The aging-regime 
temperatures for the two subsystems have become equal (thermalized case). 




Fig. 5. The susceptibility xn(M«>) and x^2{t,t w ) vs. the corresponding auto-correlation functions 
C\\(t,t w ) and C22(t,t w ) for the uncoupled, aging systems. The effective temperatures in the 
aging regime are different. 

The effective temperature indeed regulates thermalization, and we believe this 
the reason why it deserves its name. The out of equilibrium dynamics of the model 
can then be interpreted as having fast scales that are thermalized with the bath, 
and slow (aging) scales that are at a higher effective temperature. As time passes, 
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Fig. 6. The susceptibilities Xij(t,t w ) vs. the correlation functions dj(t,t w ) for the two weakly 
coupled subsystems. All the curves become parallel: also the aging regimes have thermalized. 



more modes thermalize with the bath, and the frequencies that are in the aging 
regime become lower and lower. The effective temperature shares some, but not all, 
properties with the 'Active temperature' of glass phenomenology (see Ref. 14) for 
a discussion). 
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